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Q" 1 ' Abstract 

<D ■ 

Spontaneous-fission half-lives (Ty^) of the heaviest even-even nuclei are evaluated 
and compared with their alpha-decay mode. Calculations of T-jyj are performed in 
the dynamical way with potential energy obtained by the macroscopic-microscopic 
method and the inertia tensor obtained by the cranking approximation. The alpha- 
^ ■ decay half-lives (Ty 2 ) are calculated for the same region of nuclei by use of the Viola 

\Q • and Seaborg formula. The ground-state properties such as mean square radii and 

electric quadrupole moments are also studied. 

From the analysis of T^K it is found that a peninsula of deformed metastable 
superheavy nuclei near fogSg^ * s separated from an island of the spherical super- 
heavy, around the doubly magic nucleus 298 114i84, by a trench in the vicinity of 
neutron number N=170. 
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1 INTRODUCTION 

Much progress has been made recently in the synthesis of very heavy nuclides. Deformed 
superheavy nuclides with proton numbers Z=108, 110, 111 and 112 have been discovered 
through reactions of cold-fusion at GSI |], [| [3J and hot-fusion in Dubna [f|, |5]. In 1995, 
first chemical separations of element 106 were performed |J. After these successes, the 
accurate calculations of the lifetimes of nuclei situated in the upper-end of the isotopic 
chart became a new challenge of nuclear theory. 

The objective of the present paper is to study spontaneous-fission (T-jyf ) and a-decay 
(T°/ 2 ) half-lives for even-even nuclei with proton number Z=100-114 and neutron number 
N=142-180. This relatively broad region of nuclei contains experimentally well known 
nuclides with Z<104, the deformed superheavy with Z>106 and transitional nuclei close 
to a hypothetical island of spherical superheavy elements situated around the nucleus 
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298 114 184 . For all these nuclei the action integrals describing probability of SF are min- 
imized using multi-dimensional dynamic-programming (MDP) method based on WKB 
approximation within the same deformation space. 

Much attention is paid to find the optimal deformation space. We examine a relatively 
reach collection of nuclear shape parameters {(3\, with A=2,3,4,5,6 and 8), as well as the 
pairing degrees of freedom (i.e. proton A p and neutron A n pairing gaps). The optimal 
collective space {fa, (3^ /3q, A p , A n } is found by comparison of the calculated T-jyf of Fm 
isotopes with their experimental values. 

Alpha-decay is one of the most predominant modes of decay of superheavy nuclei. 
All recently discovered superheavy elements with atomic numbers Z>107 were identified 
from their a-decay chains. Calculations of Ty 2 are easier to perform than TfE- The 
half-life for a-decay depends primarily upon the release energy Q Q , which is given only 
by appropriate difference of ground-state masses. To better characterized properties of 
the nuclei in investigated region, we also calculated their ground-state electric quadrupole 
moments and mean square radii. A part of the results of the analysis has been presented 
earlier 0, g |. 

The description of the method and details of the calculations are given in Sec. 2, the 
results and discussion are presented in Sec. 3. 

2 DESCRIPTION OF THE METHOD 

2.1 Collective variables 

Let us consider a set of A collective variables (X\, A 2 , X\) = X, then the classical 
collective Hamiltonian can be written 

1 A 

H = ~ E B kl (X)X k X t + V(X) , (1) 
1 k,i 

where V(X) is the collective potential energy. The tensor of effective mass (inertia tensor) 
Bki is symmetric and all its components may depend on collective variables. The choice 
of the collective variables is arbitrary but caution must be exercised. 

Starting from the single-particle motion of A=Z+N nucleons, a state of nucleus is 
defined by an average potential fixed by a set of parameters. These parameters are good 
candidates for collective variables. In the presented paper we used the single-particle 
Hamiltonian {H s/p ) consisted of a deformed Woods-Saxon Hamiltonian and a residual 
pairing interaction treated in the BCS approximation. In our model we discussed only 
axially symmetric deformations of H s p , i.e. a nuclear radius was expanded in terms of 
spherical harmonics Y\o(cos'd) 

^max 

fl(tf) = R (f3 x ) [1+ £ PxYm(co8&)] , (2) 

A=2 

where j3\ is the set of deformation parameters up to \ ma x multipolarity and the dependence 
of R on (3\ is determined by the volume-conservation condition. 
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Besides the deformation parameters, the other candidates for collective variables are 
parameters connected with the pairing interaction. In the usual BCS approximation we 
have two such parameters: the pairing Fermi energy (A) and the pairing energy-gap (A). 

In the presented paper we chose the proton A p and neutron A n pairing gaps as the 
additional collective variables. The significant role of these so-called pairing vibrations 
on penetration of the fission barrier has been first discussed by Moretto et al. in Ref. [[Uj 
(see, also [IT], |TJ, [TJ, |TJ). 

Finally, the set of collective variables consists of the nuclear shape parameters (f3\) 
and the pairing degrees of freedom (A p , A n ). These variables spann the multi-dimensional 
deformation space {^a} within which we shall describe a fission process. 

2.2 Inertia tensor 

The inertia tensor Bki describes the inertia of a nucleus with respect to changes of its 
shape. It also plays a role similar to the metric tensor in deformation space {X^}. Its 
components for multipole vibrations as well as pairing vibrations can be evaluated in the 
first order perturbation approximation |j| 

where |0) and \m) denote a ground state and the excited states of the nucleus with 
the corresponding energies So and £ m . For even-even nuclei the excited states can be 
identified with the two quasi-particle excitations S m = + E v . After transformation to 
the quasi-particle representation the corresponding formula takes the following compact 



form [16 1 



B H {X) = 2h 2 J2P»:(X)(E„ + E,Y l Pl v {X) , (4) 



where for the shape deformations 



and in the case of pairing degrees of freedom 

P^(A) = ^ >~^+ A ^ . (6) 

Here t> M , m m are the pairing occupation probability factors, e M are the single-particle 
energies of H s . p . and the = [(e„-A) 2 + A 2 ] 1 / 2 is the quasi-particle energy corresponding 
to state. The above expression is equivalent to a commonly used formula developed 
by Sobiczewski et al. in Ref. |T7| . 



The components of inertia tensor are strongly affected by single-particle and pairing 
effects. The relation between the energy-gap parameter A, the effective level density at 
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the Fermi energy g e ff(X) and the diagonal components of inertia tensor Bkk can be showed 
in terms of the uniform model [lITf 



Bkk ~ const .^I^-\(dH s . p JdPk)\ 2 ■ (7) 



This strong dependence of the inertia tensor on the pairing energy-gap allows to expect 
considerable reduction of the spontaneous-fission half-life values. 

2.3 Collective energy 

The collective energy V is calculated for a given nucleus by the macroscopic-microscopic 
model developed by Strutinsky [[19]: 

V = E macr ((3) + 5E sheU (f3) + 8E pail ((3, A) . (8) 

For the macroscopic part E mSuCT we used the Yukawa-plus-exponential model [pD| . The 
so-called microscopic part, consisted of the shell 5-£7 s heii and pairing SE pa i r corrections, 
was calculated on the basis of single-particle spectra of Woods-Saxon Hamiltonian [T2T 



The one-body Woods-Saxon Hamiltonian consists of the kinetic energy term T, the 
potential energy V ws } the spin-orbit term V^ s and the Coulomb potential Vcoui for 
protons: 

H ws = T + ywsft 0) + ywsff /?) + 1(1 + r 3 )V C oui(r; P) • (9) 
In the above equation 

1 + exp[dist(r; p)/a\ 

and 

KT(r; P) = -\(VV WS x p) ■ s, (11) 

where dist(f; (3) denotes the distance of a point r*from the surface of the nucleus given 
by Eq. (2) and Vo, k, a, X are adjustable constants. The Coulomb potential Vcoui is 
assumed to be that of the nuclear charge equal to (Z — l)e and uniformly distributed 
inside the nuclear surface. In our calculations we used Woods-Saxon Hamiltonian with 
the so-called "universal" set of its parameters (see Ref. [2~T|) which were adjusted to the 
single-particle levels of odd-A nuclei with A>40. 

The term 5E pauiT in Eq. (8) arises from the pairing residual interaction, which is included 
to our H sp by the BCS approximation. In the presented paper we used the pairing 
strength constants: G Z A = 13.3 + 0.217(N - Z) and G N A = 19.3 - 0.080(N - Z) for 



protons and neutrons, respectively, which are taken from Ref. [22 
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2.4 Lifetimes for alpha— decay 



For calculation of alpha-decay half-life we employ the phenomenological formula of Viola 
and Seaborg [|3[ 



(12) 



logT^ [yr] = (aZ + b)(Q a /MeV)- L/2 + (cZ + d) - 7.5 , 

where Z is the atomic number of the parent nucleus and Q Q is the energy release obtained 
from the mass excesses 



Q a (Z, N) = M{Z, N) - M(Z - 2, N - 2) - M(2, 2) . (13) 
The values of parameters: a=l. 66175, fe=-8.5166, c=-0. 20228 and <i=-33.9069 in the 



above formula were taken from Ref. |[24|| . It should be noted that, the uncertainties in 
the calculated a-decay half-lives due to their phenomenological character are far less than 
uncertainties in the calculated SF half-lives. 



2.5 Lifetimes for spontaneous— fission 

The spontaneous-fission half-life is inversely proportional to the probability of penetration 
through the barrier 



T m = — 75- (14) 



In 2 1 

Where n, in the above formula, is the number of "assaults" of the nucleus on the fission 
barrier per unit time. The number of assaults is usually equated to the frequency of zero- 
point vibration of the nucleus in the fission degree of freedom and for a vibrational fre- 
quency of hu =lMeV , assumed in this paper, n ~ 10 20,38 s _1 . Using the one-dimensional 
WKB semi-classical approximation for the penetration probability P one obtains 

1 n-28.04 

Tf/f [yr] = -j— [1 + exp 2S(L)} , (15) 

where S(L) is the action-integral calculated along a fission path L(s) in the multi- 
dimensional deformation space {X\} 

S(L) = £ 2 [^B eff (s)[V(s) - E^ds. (16) 
An effective inertia associated with the fission motion along the path L(s) is 

BM S ) = Y.B« d -^ d -^, (17) 

where Bki are the components of the inertia tensor. 

In the above equations ds denotes an element of the path length in the {X\} space. 
The integration limits Si and S2 correspond to the classical turning points, determined by 
a condition V(s) = E, where E = V(X®) + 0.5 Tiujq denotes the energy of the fissioning 
nucleus in MeV (calculated in the ground-state). 
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2.6 Calculations technique 



Dynamic calculations of Tjyf mean a quest for the least-action trajectory L m j n which 
fulfills a principle of the least-action 5[S(L)] = 0. To minimize the action-integral (|16|) 
we used the dynamic-programming method. Its application to fission was first developed 



by Baran et al. (see e.g., p5| and references cited therein). 



In contrast to the method used by Smolariczuk et al. in Ref. |26], \Z7\ , where only two 
coordinates (/?2 and fi^ have been handled dynamically and the remaining degrees of free- 
dom have been found by minimization of the potential energy V, in our multi-dimensional 
dynamic-programming (MDP) method all coordinates are treated as independent dynam- 
ical variables. 

Figure 1 demonstrates how our model works. Since the macroscopic-microscopic 
method is not analytical, it is necessary to calculate the potential energy and all compo- 
nents of the inertia tensor on a grid in the multi-dimensional space spanned by a set of 
deformation parameters {X\}. We select one coordinate Xo from this set. This coordi- 
nate (e.g. elongation parameter) is related in a linear way to the fission process. In Fig. 1 
to each point X correspond a plane representing the rest of the collective space {Xa-i}. 




{Xj=x e{x^} 

Figure 1: Diagrammatic illustration of the MDP method. In the multi-dimensional 
deformation space {X\} we select the coordinate {Xq} which is related in a linear way to 
the fission process. The points si and S2 correspond to entrance to the barrier and exit 
from the barrier, respectively. See text for details. 

To find the least-action trajectory L min between the turning point si and s 2 we proceed 
as follows. First, we calculate the action-integrals from the entrance point under the 
barrier s\ to all grid points in the nearest plane at Xq = 1. In the next step we come to 
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the plane at Xo = 2 and from each grid point in this plane calculate the action-integrals 
to all grid points in the plane at X = 1. The trajectories started from each grid point 
at Xq = 2, passing through all grid points in the plane at X = 1 and terminated in 
the point si, form a bunch of paths. From each such a bunch we choose the path with 
the minimal action-integral and bear it in mind. At the end of this step we have the 
least-action integrals along trajectories which connect the starting point si with all grid 
points in the plane at X = 2. Next, we repeat this procedure for all grid points at X = 3 
and again we obtain all the least-action-integrals along trajectories starting from point 

51 with ends at each grid point in the plane X = 3. We repeat it until we reach the 
n— th plane, the last one before the exit point from the barrier s 2 . Finally, we proceed to 
the last step of our method, where we calculate action-integrals between the exit point 

52 and all grid points situated on the last plane at Xo = n; the minimal one among them 
corresponds to the searched trajectory of the least-action-integral L min . 

If we denote a number of grid points on each Xj (i=l,2,...,A-l) axis by n^, then the 
whole number of trajectories examined in MDP method is equal to {n\ -ni • ... • n\_i) n . Up 
to now, our calculations are carried out in a maximum of four-dimensional deformation 
space in view of an enormously large computational time (and disk space) required for 
preparing (and storing) input data with potential energy and 1/2A(A + 1) components of 
symmetric inertia tensor for each of (ni ■ n 2 ■ ... ■ n\_i)-n grid points. 

Calculations are performed in various four-, three- and two-dimensional deforma- 
tion spaces spanned by selected shape parameters (/3\, with A=2,3,4,5,6 and 8) and two 
pairing degrees of freedom (A p , A n ). For /3-shape parameters we used grids with steps 
A/3 2 =A/3 3 =0.05 and A/3a=0.04 for A=4,5,6 and 8; in the case of pairing energy-gaps grids 
steps are equal 0.2 MeV. In our calculations a quadrupole deformation /3 2 plays a role of 
the coordinate X in Fig. 1. 

3 RESULTS AND DISCUSSION 

3.1 Optimal Multi— dimensional Deformation Space 

The experimental values of the spontaneous fission half-lives of nine even-even Fm iso- 
topes (N = 142, 144, 158) form approximately two sides of an acute-angled triangle 
with a vertex in N = 152. This strong nonlinear behaviour of TfE vs. neutron number 
N, due to an enhanced nuclear stability in the vicinity of deformed shell N=152, provides 
good opportunity for testing theoretical models. 

To find the proper deformation space for description of the fission process we examined 
three different effects: the effect of the higher even-multipolarity shape parameters /3q and 
/3s, the role of the reflection-asymmetry shape parameters (3% and (3$, and the influence 
of the pairing degrees of freedom A p and A n . 

The following conclusions can be drawn from our previous dynamical analysis of Tf^ 
for Fm even-even isotopes, see Ref. 05 H 0- In the case when the (3^ deformation 
parameter is added to our minimal two-dimensional space {/? 2 , (3^} we observed an increase 
of fission lifetimes by one to four orders of magnitude. The contribution of parameter (3$ 
to TfK is negligible. 
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The deformations with odd-multipolarity A=3,5 do not change SF half-lives. The 
reason of this lies in the dynamical treatment of the fission process. The parameters fa 
and /3 5 reduce the width of a static fission barrier, however the effective inertia -B e //, Eq. 
(|17|), along the corresponding static path is larger than along the dynamic one, where 
/?3 and /?5 are almost equal to zero. One can say, that the static path (corresponding to 
minimal potential energy) is "longer" than the dynamical one, for which /53=/?5=0. The 
above conclusions are in agreement with those published in Ref. [^J . 

The pairing degrees of freedom A p and A n reduce SF half-lives for Fm isotopes with 
N>152 for about 3 orders of magnitude and considerably improve theoretical predictions 
of T-jyf ■ This effect is due to strong dependence of the inertia tensor upon pairing energy- 
gap, as it was shown in Eq. (|7|). 

Finally, we can conclude that the optimal deformation space for description of the 
fission half-lives of heavy nuclei is /3 2 , /3 4 , /3 6 , A p and A n . 




Figure 2: The electric quadrupole moments of the even-even nuclei with atomic numbers 
Z= 100-1 14, plotted as a function of the neutron number. 

On account of computational limitations mentioned above, we can only perform calcu- 
lations in a maximum of four-dimensional deformation space. So, we decided to define a 
correction to SF half-lives, which arises from pairing degrees of freedom, as the difference 
between TfE calculated in four-dimensional space {/3 2 , /3 4 , A p , A n } (when pairing degrees 
of freedom are treated as dynamical variables) and the one calculated in two-dimensional 
space {fa, At} (where pairing energy-gaps are treated in the stationary way- i.e. by 
solving the BCS equations): 
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Calculations of T\ 




obtained in three-dimensional space {(3 2 , 04, P&} with the pairing correction 8T^(/S. P , A n ). 
3.2 Ground-state properties 



Figure 3: The mean square radii of the even-even superheavy nuclei, plotted as a function 
of the neutron number. 

In the present study of superheavy the even-even nuclei with atomic numbers Z=100- 
114 and neutron numbers N=142-180 are considered. First, we present the results related 
to ground-state (GS) properties. The GS properties were calculated in the equilibrium 
point found for a given nucleus by minimization of its potential energy with respect to 
P 2 , Pa and p 6 degrees of freedom. 

In Fig. 2 we plot the electric quadrupole moments calculated with following formula 



where v 2 v is the BCS occupation factor corresponding to proton single-particle state \v) 
in the equilibrium point. 

Almost all nuclei have distinct prolate deformations. And with an increase in the 
neutron number the Q 2 values show a regular decrease except at N=162-164, where a 
slight discontinuity in this behaviour can be seen for nuclei with atomic number Z>106. 



39.5 
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(19) 



9 



The mean square charge radii (MSR), for the same region of nuclei, are plotted as a 
function of neutron number in Fig. 3. For calculations of MSR we use the usual formula 



< r 2 >= - ( v | r 2 | v ) v 2 + 0.64 fm 2 

Z u =p 



(20) 



where the last term is due to finite range of proton charge distribution. 

One observes a rather regular dependence of mean square radii on both neutron and 
proton number. However, as previously, close to N=162-164 one can see local maxima 
in MSR curves, particularly for nuclei with Z>106. This means that Coulomb repulsion 
energy for these nuclei is locally smaller, then they are more stable. 

3.3 Spontaneous— fission versus alpha— decay 
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Figure 4: Contour map of the logarithm of the spontaneous-fission half-life (given in 
years) for nuclei shown in Fig. 2 calculated in {P2, Pa, Pe} deformation space and corrected 
by the effect of the pairing degrees of freedom, Eq. (18). 



Figure 4 shows the results of the spontaneous-fission half-lives calculation, according 
to the MDP method, performed in {fa, Pa, P&\ deformation space with pairing correction 
5T^f(A p , A n ), Eq. for the nuclei shown in Fig. 2. 

Two very specific effects can be observed on the contour map of T-gf plotted as a 
function of neutron and proton numbers. One can see an enhancement in the SF half-life 
values at N=162 followed by a diminution at N=170. 

The enhancement in nuclear stability near the deformed shell N=162 allows the ap- 
pearance of a peninsula of deformed metastable superheavy nuclei. The local maximum 
of the T^yf values is centered at the nucleus 2 Q6Sgi62 (2-5 h). 
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In the vicinity of neutron number N=170 one observes the opposite behaviour. Here, 
the SF half-life values form a trench. This trench separates the peninsula of the deformed 
superheavy nuclei from an island of spherical superheavy elements around the doubly 
magic nucleus 298 114 184 . The local minimum is obtained for nucleus ^No^o (10 /is). 

We found also that the T®E values of two heaviest nuclei (considered in the presented 
paper) 292 114 178 and 294 114 180 are comparable with those of the most stable Fm isotopes 
with neutron numbers N=152 and 154 (~ 100 yr). 

In Fig. 5 we plot the alpha-decay half-lives (given in years) estimated by use of the 
Viola-Seaborg relationship with set of constants from Ref. [P^| . The contour plot of 
the Ty 2 forms a relatively regular surface descending steeply in the direction where the 
proton number tends to increase and the neutron number tends to decrease (upper-left- 
hand corner in the plot). 



114 
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Figure 5: Contour plot of the logarithm of the alpha-decay half-life Ty 2 (given in years) 
obtained from the Viola-Seaborg systematics. The alpha-decay energy was calculated 
with {j3\}, A = 2, 4, 6 deformation parameters. 



The surface of the Ty 2 values shows an evident protuberance at N=162, which demon- 
strates unambiguously the magicity of this neutron number. The shell effect at N=152 is 
very weak and disappears practically for nuclei with atomic number Z>104. 

The results presented in Fig. 4 and 5 as well as conclusions drawn from them are 
generally similar to those recently published by other groups employing the macroscopic- 
microscopic method, Ref. |27], |2S|, ^ [HJ. 

To compare the SF and a-decay modes we show in Fig. 6 the logarithm of the total 
half-life T^ +a resulting from both these modes. If we examine contour maps with T^ +a 
and TfE we can notice some minor differences but the global behaviour of both these 
quantity stays unchanged. 
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The dark shadowed areas in Fig. 6 show the regions of nuclei in which the a-decay 
mode predominates. The light shadowed area corresponds to the intermediate region of 
nuclei, where the probabilities of the SF and a-decay processes are approximately equal 
(i.e. the region where the values of Ty 2 and differ up to one order of magnitude). The 
black solid curve inside this area connects nuclei for which probabilities of both considered 
modes are the same. 

Thus, one can observe that the region of increased a-decay activity separates two 
areas with predominating SF activity in a diagonal manner, in Fig. 6. It is worth noting 
that the upper area is almost inaccessible due to extremely short lifetimes (< 1 /j,s). 



logT^ +cx [yr] 



114 




100 
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Neutron Number N 

Figure 6: Logarithm of the total half-life T x K^ 01 resulting from SF (Fig. 4) and a-decay 
(Fig. 5) modes. The dark shadowed areas show the regions of nuclei, where the a-decay 
mode is a dominant one. The light shadowed area corresponds to the intermediate region, 
where the probabilities of SF and a-decay modes are approximately equal. 

The total half-life values (equal to Ty 2 ) for two heaviest nuclei, considered in this 
paper, 292 114i7g and 294 114igo are found larger than 1 yr. This is in agreement with 
results obtained recently in the fully selfconsistent microscopic nonrelativistic Hartrre- 
Fock-Bogoliubov approach, Ref. [plj. 
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